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USING RESIN AGE FACTOR TO OBTAIN 
MEASUREMENTS OF IMPROVED ACCURACY OF ONE OR MORE POLYMER PROPERTIES 

WITH AN ON-LINE NMR SYSTEM 



Cross-Reference to Relat ed Patents and Applications 

This application is a continuation-in-part of pending U.S. patent application serial number 
08/586,559 filed on January 16, 1996 which is a continuation-in-part of pending U.S. patent 
application serial number 08/491,632 filed on June 19, 1995 which is a continuation-in-part of the 
5 following two U.S. patent applications: (1) U.S. patent application serial number 08/370,862 filed 
on January 10, 1995 which is a continuation-in-part of U.S. patent application serial number 
08/226,06 1 filed on April 1 1 , 1 994 (now abandoned) which is a continuation of U.S. Patent 
No. 5,302,896 issued to Dechene et al. on April 12, 1 994; and (2) U.S. patent application serial 
number 08/371,091 filed on January 10, 1995 which is a continuation-in-part of U.S. patent 
1 0 application serial number 08/226,024 filed on April 1 1 , 1 994 (now abandoned) which is a 
continuation of U.S. Patent No. 5,302,897 issued to Dechene et al. on April 12, 1994. 

Also, this application is related to the following U.S. patents: 5,015,954 issued to Dechene 
etal. onM.ay 14, 1991; 5,049,819 issued to Dechene et al. on September 17, 1991; 5,319,308 
issued to Dechene et al. on June 7, 1994; 5,367,260 issued to Dechene et al, on November 22, 
1 5 1 994; and 5,408, 1 8 1 issued to Dechene et al. on April 1 8, 1 995 . 

Each of these U.S. patents and patent applications is of common assignment with this 
application, and the disclosures of all of the above-listed patents and patent applications are 
incorporated herein by reference. 
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Field of the Invention 

The present invention involves using pulsed nuclear magnetic resonance (NMR) techniques 
to measure physical properties (e.g., xylene solubles, density, rubber/oil content, and flow rates 
such as melt index, flow rate ratio, melt flow) of polymer materials (e.g., rubber and plastics such 
as polypropylene, polyethylene, and polystyrene) in a real-time, on-line Environment. 

Background of the Inwnt mp 
The pulsed NMR techniques described herein, and in the above-identified related patents 
and applications, involve the use of a pulsed burst or pulse of energy designed, to excite the nuclei 
of a particular nuclear species of a sample being measured. The protons, or the like, of the sample 
are first precessed in an essentially static magnetic field. The precession thus ; is. modified by the 
pulse. After the application of the pulse, there occurs a free induction decay (FID) of the 
magnetization associated with the excited nuclei. That is, the transverse magnetization associated 
with the excited nuclei relaxes back to its equilibrium value of zero, and this relaxation produces a 
changing magnetic field which is measured in adjacent pickup coils. A representation of this 
relaxation is the FID waveform or curve. 

An NMR system described herein, and in the above-identified related patents and 
applications, takes measurements of physical characteristics of polymer materials (e^g , rubber, 
plastics, etc.) and relates those measurements back to, for example, flow rates (e.g., melt index), 
crystallinity, composition, density, and tacticity by performing the following analysis methods. 
The NMR system is first calibrated with known samples by determining the physical types, . 
properties, and quantities of the target nuclei in each known sample. Unknown samples are then 
introduced into the NMR system, and the system determines the physical types, properties, and 
quantities of the target nuclei in each unknown sample based on the calibration information. 

The analysis methods performed by the NMR system involve decomposing an FID curve 
associated with a known sample into a sum of separate time function equations. Useful time 
function equations include Gaussians, exponentials, Abragams (defined herein as Gaussian 
multiplied by the quantity sin(cot) divided by cot), modified exponential (defined herein as Ce z 
where C is a constant, z=(kt) a , and a lies between 0 and 1 or 1 and 2), modified Gaussian (defined 
herein as Gaussian multiplied by the cosine of the square root of cot), and trigonometric. 

The coefficients of the time function equations are derived from the FID by use of a 
Marquardt-Levenberg (M-L) iterative approximation that minimizes the Chi-Squared function. 
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This iterative technique is well-known in the art. Other known-in-the-art iterative techniques can 
be used instead ofM-L including a Gauss-Jordan technique, a Newton-Raphson technique, and/or 
a "steepest descent" technique. 

From the time functions, a set of parameters is calculated. Some of these parameters are 
• ratios of the y-axis intercepts, decay times (or time constants of decay) for each of the time 
functions, and the cross products and reciprocals thereof. The sample temperature may form the 
basis for another parameter. 

: Statistical modeling techniques are then used to select a subset of these parameters to form 
a regression equation or model. Regression coefficients are then computed for this parameter 
subset. These regression coefficients represent the regression model which relates the parameter 
subset to the types, properties, and quantities of the target nuclei in the known sample. 

After the NMR system has been calibrated with one or more known samples, unknown 
samples can be introduced thereinto. 

When an unknown sample 1'is int^duced into the calibrated system, the FID 
associated with the uhknbwn.anipie is decomposed into a sum of separate time function 
equations. The coefficients of the time function equations are derived from the FID by use of the 
iterative M-L technique: From the time functions, parameters are calculated. The parameters are 
then "regressed" Via the regression model to yield the types, properties, and quantities of the 
target nude, in the unknown sample. That is/ the measured parameters from the FID of the 
unknown sample are used with the regression model, and the types, properties, and quantities of 
the target nuclei in the unknown sample are determined. This informauon is related to a property 

Otte s^l^^m *^, xylene solubles, or flow rates (e.g., melt index, flow rate 
ratio, and/or melt flow). 

The regression model typically is a muiti-dimensional regression equation or model and it 
may be hnear or non-linear. Because it is muiti-dimensional, it may not be graphically 
represented. As an example of a regression technique, consider that the grade point average of 
each of the students at a college were related to that student's SAT score and high school 
standmg; forming a three dimensional space. The line formed is a "regression line" which may be 
graphed. A new student's grade point average may be predicted by placing the new student's 
SAT score and high school standing on the "regression line" and "reading" the grade point 
average. 
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Melt index (MI) has been defined for polyethylene as the flow rate obtained under 
condition 190/2. 16 (note 19, ASTM, or American Society for Testing and Materials, No. D1238- 
90b). In general, MI is a measure of viscosity. Flow rate ratio (FRR) has been defined for 
polyethylene as a dimensionless number derived by dividing the flow rate at condition 190/10 by 
the flow rate at condition 190/2.16 (paragraph 8.3, page 396, ASTM Ne. D1238-90b)! In some 
cases, the logarithm of FRR is used instead of FRR. Melt flow (MF) has been defined for 
polypropylene as the flow rate at condition 230/2. 16 (ASTM No. D1238-90b). 

The age of a sample of polymer material ("resin") under test may have an impact on the 
accuracy of the property measurement made with the NMR machine. However, there generally 
have not been any attempts to understand or account for "resin age" in measurements made with 
an real-time, on-line NMR machine. . ; ' 
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Summarv of the Invention 
It has been found that the age of a sample of resin impacts the accuracy of measurements 
made of polymer properties (e.g., MI or average molecular weight) with a real-time, on-Jine 
nuclear magnetic resonance (NMR) machine. The term "age" or "resin age" as used herein is the 
time interval between the actual format.on of the polymer and the actual NMR measurement of a 
sample of the polymer. 

It is an object of this invention to improve the accuracy of the estimations of polymer 
properties (e.g., MI) produced by an on-line NMR system for various polymer materials. 

It is also an object of this invention to use a "resin age" factor to enhance the predictive 
accuracy of the.NMR system and thus obtain more accurate pro P ert y -of-interest (e.g., xylene 
solubles, density, rubber/oil content, or flow rates such as MI, FRR, and MF) measurements of a 
Plastic (e.g., polypropylene, polyethylene, or polystyrene) or any other polymer including rubber. 

The invention relates to "on-line" techniques and systems that do not require any external 
input, other than the unknown sample itself, in order to determine accurately the sample's 
property of interest. 

In accordance with the invention, process data for generating system models of enhanced 
predictive accuracy includes a resin age factor which accounts for time dependent aging 
phenomena and resultant changes in values of polymer crystallinity, melt index, density, and other 
polymer properties of interest. The resin age factor typically is a function of one or more process 
parameters such as elapsed cycle time prior to sample measurement, temperature of the sample at 
tune of measurement, and sample form. By using a resin age factor, the model can be simplified 
by reducmg the number of discrete input terms required and overall predictive accuracy can be 
enhanced. 

The foregoing and other objects, aspects, features, and advantages of the invention will 
become more apparent from the following description and from the claims. 
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Brief Description of the Drawing s 

In the drawings, like reference characters generally refer to the same parts throughout the 
different views. Also, the drawings are not necessarily to scale, emphasis instead generally being 
placed upon illustrating the principles of the invention 

FIGS. 1A, IB and 1C are block/schematic diagrams of a prefen^d embodiment of a pulsed 
NMR system according to the invention which is suitable for measuring various properties of a 
range of polymer materials. ' ' 

FIGS. 2A-2D are diagrams of other embodiments for handling samples. 

FIG. 3 is a diagram showing a sensor for measuring the sample temperature. 

FIG. 4 is a graphical representation of a free induction decay (FID) curve and its 
component curves measured for polyethylene. " 

FIG. 5 A is a flowchart of steps performed according to the invention to determine 
properties of a polymer sample. J v; * 

FIGS. 5B and 5C are graphs of NMR measurements bf density versus NMR 
measurements of melt index (MI) showing various plastic grides (polyethylene in this case). 

FIGS. 6 A and 6B are flowcharts of steps, more detailed than those shown in FIG. 5 A, 
performed in accordance with the invention to determine properties of an unknown polymer 
sample. . ' : •••• 

FIG. 7 is a flowchart of the steps to establish an effective industrial measurement . 

FIG. 8 is a flowchart of steps performed in accordance With an alternate embodiment of 
the invention to determine properties of an unknown polymer sample. 

FIG. 9 is a schematic depiction of an exemplary riieural network useful for generating 
predictive models in accordance with another alternate embodiment of the invention. 

FIG. 1 0 is a flowchart of steps performed in accordance with the neural network derived 
model depicted in FIG. 9 to determine properties of an unknown polymer sample. 

FIG. 1 1 is a flowchart of steps performed in accordance with the invention utilizing an 
exemplary embodiment of a resin age factor cycle time parameter. 

FIG. 12 is a flowchart of steps performed in accordance with an embodiment of the 
invention employing a resin age factor to determine properties of an unknown polymer sample. 
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Description 

FIG. 1A shows transverse and cross sections, with block diagram inserts, of an NMR 
apparatus and method where the present invention may be used to advantage. An industrial 
process line D>L has material flowing as indicated by arrow A. Some of the material is captured 
by a probe P and fed through an inlet line LI to a sample region S 1 . .The region is defined by a 
tube 98 typically about 30 cm long made of an essentially non-magnetic, nonconducting material 
which does not itself generate substantially interfering free induction decay (FID) signals. The 
tube material can be, for example, glass, certain ceramics, certain plastics, or hybrids thereof. The 
sample region is defined between inlet and outlet valves VI and V2. Gas jets J are also provided. 
These jets are pulsed pn/ofF repeatedly to agitate fluent sample materials during sample admission 
and expulsion. A region S2 is the critical portion of the sample. It is surrounded by a sample coil 
100 tuned to resonance and driven by a tuning circuit 102. and related transmitter/receiver 
controller 104. Grounded loops 101 are Lenz Law shields which are provided above and below 
coil 100 to help shape the field .of coil l Oq. j.e, to contain the field established by an excitation 
pulse. The controller 1 04 includes an on-board microprocessor and required power supply 
elements, memory, program and I/O (input/output) decoding suitable to interconnect to the 
hardware shown and to an external microcomputer 1 06 with associated memory, a keyboard 108, 
a monitor (or other display) 1 10, a recorder 112, and/or a process controUer-1 14 (to control the ' 
process at IPL). The operator initiates and pontrols operation from the keyboard 1 08 and the 
resulting data and signals are subsequently shown on the display 1 10 and utilized in the recorder 
1 12 and/or the controller 114 The computer 106 also controls instrument operation conditions. 

The region S2 of tube 98 and coil 1.00 are in a static, but adjustable, crossing magnetic 
field defined by a magnetic assembly 1 16 which comprises a yoke 1 18, pole pieces 120, 
surrounding Helmholtz coils 124, and a coil current generator 1 17. The critical sample' region S2 
of the tube 98 and magnet are contained in a metallic (but non-ferromagnetic) box 126 with highly 
thermally conductive face-plates 128..and internal partitions 130 and overall mass related to each 
other to minimize harmonics and other interferences, with a signal emitted from coil 1 00 to a 
sample and/or returned from the sample for pick-up by coil 1 00 and its tuned circuit i.02 and 
transmit/receive controller 104. 

The magnetic assembly, including yoke 118, and other parts therein as shown on FIGS. 
1A and 1C, is in turn contained in an environmental control chamber 132 with optional inert gas 
fill and purge controls (not shown), an internal electrical heater 134, a motor M driving fan 136 
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and a temperature sensor 138 in the air stream whose temperature is reflective of the temperature 
at pole pieces 120 and in the region surrounding the air curtain containment jacket 137. 
Additionally, there is a heated chamber shown in FIG. IB where the sample is heated prior to 
introduction into the NMR sample chamber. A thermal controller 140 processes temperature 
signals from the sensor 138 to adjust heating/circulation at the heater/fay 134/136 as a coarse 
control and to adjust current through the Helmholtz coils 1 24 at magnet pole pieces 1 20 as a 
sensitive and fast fine control, as well as implementing general control instructions of computer 
106. Further thermal stabilization is provided by a temperature controlled air curtain consisting of 
temperature controlled air 135 circulating inside an air containment jacket 137, completely 
surrounding the sample region under NMR measurement. The temperature of the air curtain is 
controlled by the heater thermal control 1 3 IB and temperature sensor 13 1 A via 1 3 1 D according 
to control parameters received from CPU 105-2 (I/O). Sample temperature is measured by 
temperature sensor 133 from which an external computer system 106 determines the desired air 
curtain temperature set point to be conveyed to heater thermal control 13 IB. ''' 

While the terms "air" and/or ! 'air curtain'' are used herein, it should be understood that 
other gas or liquid environments can be utilized. 

FIG. IB is a schematic outline of the heating chamber 1 70. The sample is fed (as indicated 
at arrow 172) into the chamber. There is a heating element 1 74 surrounding the chamber with its 
temperature set about 10 or so degrees centigrade higher than a mobility enhancing temperature 
known for the specific material. The temperature may be looked up in the below-mentioned text. 
The mobility enhancing temperature of a polymer material (e.g. , a plastic such as polypropylene, 
polyethylene, and polystyrene) is defined herein as follows. It is the temperature that enhances 
precision and reliability for property measurements such as measurements of viscosity, molecular 
weight, melt index (MI), flow rate ratio (FRR), and melt flow (MF). Molecular weight, MI, FRR, 
and MF are closely related. Polymers generally have higher mobilities and an enhanced NMR 
response when at the mobility enhancing temperature. The mobility enhancing temperature is a 
minimum or threshold temperature in that there is a temperature range extending above the 
mobility enhancing temperature where the benefits occur without the sample handling problems 
associated with molten polymers. The mobility enhancing temperature is further defined and 
explained in the next paragraph. 

The mobility enhancing temperature is the temperature at or above the "glass transition" 
temperature (T 8 ) for amorphous polymers (e.g., polystyrene), and at or above the "crystalline 
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transition" temperature (T 0 ) but below the "melting temperature- (T») for both crystalline and 
semi-crystalline polymers. T g> T„, and T m are defined herein as temperatures where the following 
physical phase transitions take place. At the glass transition temperature (T s ), the fractional free 
volume for amorphous polymers increases appreciably and the polymer chains start undergoing 
rapid, almost isotropic motion. Amorphous polymer at this temperature has not melted and still 
retains a solid-like appearance to the lay observer. This state is sometimes referred to, in the art, 
as "rubbery". As a result, the NMR parameters, including time constants (T2's), are influenced. 
This enables better correlation with, for example, viscosity and MI or MF. Forboth crystalline 
and semi-crystalline polymers, it is.necessary to carry out NMR measurements at or above the 
crystalline transition, temperature (T u ) in order to influence the NMR parameters and obtain better 
estimations of, for example, viscosity and MI or MF In the case of a semi-crystalline polymer, 
the melting and crystalline transition temperatures are greater than the amorphous glass transition 
^^ a ^V::^W8.^^P«?^efnperature to such temperatures as just described enhances 
NMR precision and reliability for at least the measurements of the above-identified properties and 
generally for other, phys^ ; property measurements performed via NMR techniques. The specific 
temperatures as they relate to specific polymers are found in various handbooks and textbooks 
SUGh 35 Textbook ofPoly mer Science by Fred W Billmeyer, Jr. (Wiley-Interscience, a Division of 
John Wiley and $pns, Inc., 2nd ed. s 1 971). Again, these temperatures generally are collectively 
referred to herein as the mobility enhancing temperature The actual temperatures used for 
various types of polymers are determined experimentally. 

The polymers of particular interest herein are polypropylene, polyethylene, and 
polystyrene. In general, two types of polystyrene exist: "crystal" and "extended". Extended 
polystyrene is also referred to as "impact' 1 or -rubber extended". The property of pplystyrene, 
regardless whether it is crystal or extended, of most importance herein is MI. Crystal polystyrene 
typically has oil added to it, and thus the percentage of oil (oil %), in addition to MI, can be a 
property of interest to be measured when the sample is a crystal polystyrene. Extended 
polystyrene typically also has oil and/or rubber added, and therefore MI, oil %, and the percentage 
of rubber (rubber %) are measurements of interest for extended polystyrene samples. Other 
polymers that can be tested accurately with the NMR system according to the invention include 
Ethylene Propylene Rubber, Acryonitrile Butadiene Styrene (ABS) plastic, and Ethylene Vinyl 
Acetate (EVA) copolymer. 
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Still referring to FIG. IB, the sample fed (as indicated at arrow 172) into the chamber 
remains there long enough to ensure that the sample has achieved the proper temperature. The 
temperature need only be within a few degrees of an actual set point. To further ensure proper 
and uniform heating, a heated temperature controlled stream of nitrogen gas 176 or air is 
introduced via a port 178 in order to fluidize the material. After two (o^more) minutes the 
sample is introduced into the NMR measuring chamber S2 by opening the valve VI . The gas is 
introduced in a direction and with a pressure to create a turbulence in the material that induces 
movement of particles and homogenization of material temperature. 

Other embodiments for handling (e.g., heating, measuring, and discarding) the samples are 
possible. For example, referring to FIGS. 2A-2D, the sample is heated inside the probe region 
prior to the measurement. In the embodiment illustrated in FIG. 2A, the sample is contained in a 
sample holder or "basket" 300. The holder 300 itself may be positioned on an air-cylinder driven 
mechanical sliding mechanism 302 for physically removing the sample after each measurement by 
lowering the basket 300 and removing the sample mechanically and by air jets. This arrangement 
is particularly suitable for those samples that stick together during the heating time and cannot be 
removed by simple air pressure and gravity. In another embodiment illustrated in FIG, 2B, a 
sample holder 304 is permanently fixed in position in the center of tne probe and a valve 306 is 
located at the bottom of the sample holder 304. The sample is removed by lowering the valve, 
and the sample is forced out by the heating air driven from the top. The controlled heating air 
enters the actual probe chamber from the top, and it is forced around and through the sample for 
the case of pelletized material (FIG. 2C) or around and through a central vent pipe for the case of 
powdered material (FIG. 2D). The actual heating time is pre-determined in each case for 
adequate heating to the desired temperature. Sample temperature is determined by a 
thermocouple located just above the sample 

In FIG. 1C, all the electronics and computers and other elements of the system controller 
are grouped for simplicity of description into one block 106-1 . The temperature of the sample 
chamber is measured by a temperature sensor 13 3, such as a silicon junction device or 
thermocouple, embedded in the chamber wall. FIG. 3 shows another preferred embodiment 
where the sensor is an infrared (heat) sensor with a window through the sample chamber wall 
whereby heat radiation impinges the sensor. The sensor signal is fed through an I/O interface to 
the computer portion of the controller 106-1 . In another preferred embodiment, the I/O interface 
may contain an analog-to-digital (A-to-D or A/D) converter separated from the computer and/or 
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input to the computer via some other communication port, e.g., a serial port. The computer can 
interrogate the sensor and measure the temperature at any time, typically at the start of the NMR 
measurement time period and again at the end. In response to a temperature change during the 
NMR measurement, the system adjusts the temperature through the heater 13 IE. The 
temperature of the heated air 135 is measured by the sensor 131 A via 131D which is placed in the 
air path 135. The controller 106-1 maintains the sample temperature' constant through the heated 
air within the air curtain. These two means give the entire system the ability and capacity to 

maintain the sample temperature essentially constant while the environment air around the NMR 
system changes. 

As discussed before, to enable better correlation of the NMR data with viscosity and MI 
or MF (average molecular .weight), amorphous polymers are measured at or above the glass 
transition temperature, and crystalline, and .semi.crystalline polymers are measured at or above the 
crystalline transition tempore .but ^low the melting temperature (i,e ; , polymer samples are 
measured at the mobility enhancing temperature as defined above). 

Referring b^k ^ nG 1 A, the strength, consistency, and constancy of the magnetic field 
between poles 120 in. the region S2 pfthe sample is thus controlled by a uniform base magnetic 
field modified by a small correction generated by Helmholtz. coils 124 in the entire region S2. The 
Helmholtz coils 124 are energized by the coil current controller 1 17 to accurately trim the final 
magnitude of the field in which the sampje Js placed. This field is the vector addition of the fields 
due to the magnet poles 120 and the Helmholtz coils 1 24. . The controller 1 1 7 sets the current 
through the Helmholtz coils 124 using current generators. Thecoils 1 24 are wound around the 
magnet pole pieces such that the magnetic field created by the current in the coils 124 can add to 
or subtract from the field created by the magnet pole pieces. The magnitude of the current 
through the coils 124 determines the strength of the field added to or subtracted from the field due 
to the magnet pole pieces (and related yoke structure) alone. 

The actual determination of the current through the Helmholtz coils is accomplished by 
carrying out the magnetic energy and resonance, techniques hereinafter described in preliminary 
runs and adjusting Helmholtz current until the maximum sensitive resonance is achieved. 
Operating the system this way at resonance or on resonance is the preferred mode of system 
operation, in another embodiment, the Helmholtz current is set to operate the system off 
resonance by a given offset of about 0. 1 -3 KHz. 
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The major elements of electrical controls are tuner 102, including coils 100 and 101 and 
variable capacitors 102-1 and 102-2, resistor 102-3 and diodes 102-4 and constructed for tuning 
to Q of twenty to sixty to achieve coil 100 resonance, and control 104 including a transmit/receive 
switch 104-1, a transmitter 104-2 and receiver 104-3, a crystal oscillator 104-4, gated pulse 
5 generator (PPG) 104-5, and phase shifter 104-6. The crystal provides ^nominal twenty 

Megahertz carrier which is phase modulated or demodulated by the MOD, DEMOD elements of 
transmitter 104-2 and receiver 104-3. The receiver includes variable gain amplifier elements 104- 
3 1 and 104-32 for operation. The analog signals received are fed to a high speed at least 12 bit 
flash A/D converter 105-1 and internal (to the instrument) CPU element 105-2, which provides 
1 0 data to an external computer 1 06 which has a keyboard 1 08, monitor 1 1 0, modem 1 09, recording 
elements 1 1 2, and process controller elements 114, e.g., for control of valves VI, V2 via valve 
controls 115 and/or to coil current controls 1 17, all via digital-to-analog (D-to-A or D/A) 
converters (not shown). 

The analog signal FID curve is conditioned by a Etessel filter which acts as a pre-filter and 
1 5 an anti-aliasing filter as the subsequent sampling (i.e., digitizing) is usually done at 1 0 MHz. After 
digitization, the signal may be time smoothed by a Fast Fourier transform, Savitsky-Golay, or 
other filter or filtering procedure. The combination of these filters produces a relative 
improvement in signal-to-noise ratio which generally enhances the accuracy of the NMR system. 
The excitation of coil 100 and excitation-precession of the sample's proton content and 
20 subsequent relaxation/decay produces a received FM signal that! after demodulation, controlled 
gain amplification, and A/D conversion produces the free induction decay (FID) curve. 

As indicated earlier, preferred embodiments are found in the NMR instrument whether 
operated in resonance or out of resonance, including those preferred embodiments found in the 
above-listed related patents and patent applications which have been incorporated by reference 
25 hereinto. The following preferred embodiment describes a system and method for the analysis of 
polymers, in particular, polypropylene, polyethylene, and polystyrene. 

Referring to FIG. 4, the digitized FID curve data for the material under test (e.g., 
polyethylene, polypropylene, or polystyrene) are transmitted to the external computer 106 where a 
computer program finds the best coefficients of the component curves to fit each component of 
30 the digitized FID curve. In the disclosed embodiment, there are three components: an Abragam, a 
Gaussian, and an exponential. Other embodiments have more or less than three component curves 
and other curve types such as modified exponential and/or trigonometric. The determination of 
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the types of curves which make up the FID curve is important because, once the curves are 
known, they can be extended back to a time origin (shown as AO, BO, and EO at tO) which is close 
to the center of the transmitted burst signal. This is important since there are saturation effects of 
the instrument's electronic gear which occur during and immediately after the excitation burst 
signal. During this time, measurements cannot be accurately taken, yet the region of interest 
under the curve, which is a measure of the number of nuclei in. the sample, extends from the 
immediate end of the excitation burst tq where the curve is too small to be digitized or is in the 



noise. 



Again, the entire FID curve is decomposed into component curves. Coefficients defining 
the equations of the component curves are derived by using an iterative process based upon the 
Marquardt-Levenberg (M-L) approximation technique. The derivation is applied automatically 
through a structured realization in software. This M-L technique is used to determine the 
magnitude of all the parameters, constants, frequencies, etc. which best fit the FID curve. M-L is 
an iterative technique where the entire -curve, is determined at once. The M-L iteration process 
performs the curve fitting by attempting to minimize the Chi-Squared error function (the sum of 
the squared differences between the measured data points and the data points from the derived 
equation). The results of the M-L approximation are accepted if the Chi-Squared error is small 
enough and if the number of iterations to reach an acceptable Chi-Squared error does not exceed 
a preset limit, . usually about 30. If the error is not small, enough, the M-L fitting procedure may be 
re-applied with a different set of sorting assumptions. If this re-attempt also fails, the sample is 
discarded and a new sample obtained, The M-L technique is documented in the following 
references which are herein incorporated, by reference: Ind. Annl Math v^i , 1)Pp . 431-444 by 
D.W. Marquardt, 1 963: Data Reduction ,nH a w ^, h for the phvsiral ^ rmm ^ 
McGraw Hill) Chapter 1 1 by Philip R. Bevington, 1969; and The State nf tht> An in MiimwiMi * 
Analysis (London: Academic Press, David A.H. Jacobs, 1977), chapter III.2 by J.E. Dennis. As 
applied to the measurement regime of interest herein, in a preferred embodiment of the present 
invention, the selected parameters taken from the derived curves are the y-axis intercept ratios, 
time constants, frequency terms, and other parameters described below. 

Other known-in-the-art iterative techniques which may be applied instead of, or with, the 
Marquardt-Levenberg include: Gauss-Jordan and "steepest descent" (found in the above J.E. 
Dennis reference); Newton-Raphson (known in the art); "partial least squares"; or similar 
techniques including combinations of these techniques. 
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One of the major difficulties in making use of iterative curve fitting techniques (such as 
Marquardt-Levenberg) is their tendency to reach incorrect solutions. Such solutions frequently 
(but not always) contain parameters which would imply a negative quantity of protons or an 
exponential "decay" which grows with time. These incorrect solutions lead to serious errors in 
the result found for a physical sample, for example, the density or flow properties (e.g., Ml) in 
polyethylene or the extent of tacticity or MF in polypropylene. 

Once the equation of the FID curve is known, each component curve can be extrapolated 
back to the midpoint of the excitation signal to establish the intercept of each said component 
curve. 

The resulting data utilized in the computer 106 (FIGS. 1 A and 1C) is the equation for the 
FID curve as composed of a number of component curves. Each of these curves (and their 
intercepts) has been experimentally and theoretically related to particular nuclei of interest. In 
particular, when the FID curve equation is determined, the ratios of the y-axis intercepts, the cross 
product and squares of these ratios and the decay times for each of the curve components, and the 
product temperature form one or more multidimensional models for the property. 

Calibration of the system is accomplished by first measuring a number of known samples 
and using the M-L technique to derive the model equation constants associated with each known 
sample. Various non-linear transforms may then be applied to these constants, usually with the 
goal of linearizing their relationship to the dependent (i.e., predicted) parameter. Useful non- 
linear functions include exponential, logarithmic, 1 powers, and cross products of the independent 
(i.e., measured) parameters. 

FIG. 5 A is a high-level flowchart of steps performed to determine properties of a polymer 
sample (e.g., a sample of plastic such as polypropylene, polyethylene, or polystyrene). The NMR 
system develops one or more regression equations or models for a given property (e.g., two 
models for MI, one model for density, etc.) during a calibration procedure using known samples 
(steps 1 50-1 56). In some embodiments and for some polymer properties, two or more 
regression models are developed during calibration with known samples. In one possible 
embodiment, at least two property (e.g., Ml) regression models (e.g., MIi and M2 2 ) are developed 
(step 156). A regression model for each of one or more discrete (Lei., two-valued) variables also 
are developed during calibration (step 156), and these models allow a prediction to be made about 
which of the plurality of MI regression models should be used for any particular unknown sample 
(step 164). If there are two regression models for the particular polymer property of interest, only 
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a single discrete variable regression model is needed, and the two possible values of the variable 
correspond to the two property models. That is, in the case of two property models, the single 
discrete variable model yields a value (one of the two possible values) which corresponds to the 
property model that will produce the most accurate measurement of that property for an unknown 
sample. However, if there are N regression models for the property of interest and N is greater 
than two, N discrete variable regression -models are used and each one corresponds to a different 
one of the N property models. In this N>2 case, the property model that corresponds to the 
discrete variable regression model yielding the largest value is selected as the property model that 
will produce the most accurate measurement of the property of interest for the unknown sample. 
In either case, after having selected one of the plurality of property regression models by using the 
discrete variable recession model(s) (step . 164), the selected property model is applied to data 
derived from the unknown sample to determine or estimate a value of the property of interest 
(e.g., MI) of the unknown sample (step 166). This property measurement can yield acceptable 
accuracy. In one embodiment, an acceptable MI measurement made by the NMR system has a 
standard deviationjess than 10% and typically about 5-7% or less. 

This standard deviation is the standard deviation between the MI estimation obtained by 
the NMR system and the "actual" MI, value (e.g., the MI value obtained by a manual laboratory 
ASTM method). The "actual" MI value, can be determined "off-line" in the laboratory with a 
laboratory instrument. 

..... The terete variable mentioned in the paragraph preceding the above paragraph is created 
ahead of time and prior to model making based on, for example, laboratory density a n d MI data 
for polyethylene or laboratory solubjes/tacticity and MF for polypropylene. This data can be 
gathered from various sources such as users of NMR systems which conform to the descriptions 
in the above-listed related, patents and patent applications. For example, the discrete variable can 
have a certain value (e.g., 1) for polyethylene. with,a density < 0.9225 grams/milliliter and an MI < 
1.25 and another value (e.g., 0) for any other polyethylene grade. It should be understood that 
more complex conditions might be applied; such as density and/or MI within a certain range (e.g., 
0.91 <density<0.9225 and 0.2<MI<1 .25).. This discrete variable (or variables) is referred to 
sometimes herein as INBOX because its value indicates whether or not the material is inside of the 
"box" denned by the above-identified density/Ml range. Using such a discrete variable or 
variables (i : e , INBOX,, INBOX,, ...) can result in improved polymer property measurements 
(e.g ;i ^1 measurements with a standard deviation less than 10% and typically about 5-7% or less) 
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regardless whether an unknown sample is outside or inside of a particular range or region (e.g., 
the density/Ml range identified above). That is, material for which the NMR system might 
previously have produced unacceptable property (e.g., MI) estimations (e.g., plastics falling within 
the above-identified range) may now be tested with the NMR system and an acceptable, useful 
MI estimation is achieved. 

It is important to note that the one or more "boxes" corresponding to the one or more 
INBOX variables need not be square or rectangular. In general, each "box" can have any shape, 
and each identifies a particular subset of the data Each "box" is in general a range or region. 

Referring to FIG. 5B, it was discovered that some plastic materials (e.g., polyethylene 
without hydrogen added during polymerization as indicated by circles 200, 202, and 204) fall 
within a box defined by the above-given density/Ml range, whereas other plastic materials (e.g., 
polyethylene with hydrogen added as indicated by circles 206, 208, 210, and 212) fall outside of 
the box. For the products outside of the box, the NMR system produced acceptable MI 
estimations (e.g., a standard deviation of 8:5%) by prior data analysis schemes. For the products 
inside of the box, however, the NMR system produced unacceptable MI estimations (e.g. , a 
standard deviation of about 10-12% or higher). This problem can be addressed by modifying the 
NMR system such that it now can more accurately estimate Ml, and a variety of other polymer 
properties, regardless whether the material being tested is within or outside of the box The NMR 
system automatically produces polymer property (e.g., MI) estimations which are acceptably 
accurate (e.g., MI measurements less than 10% and typically 5-7% or less) for materials falling 
both within and outside of the box. ' ' 

In general, the "grouping" or association of certain polymers to a particular model is made 
based on or determined by significant chemical properties or process conditions for that polymer 
group. 

In one embodiment, if the sample is a polyethylene sample falling within the range, a 
certain MI regression model is used. If the sample falls outside of the range, a different MI 
regression model is used. Regardless of the particular plastic being sampled, the NMR system 
automatically selects the MI model which results in the most accurate MI measurement. 

Referring to FIG. 5C, it was discovered that some plastic materials (e.g., material 
indicated by circles 220, 222, and 224) fall within a first box defined by the above-given 
density/Ml range, whereas other plastic materials (e.g., material indicated by circles 226 and 228) 
fall outside of the first box but inside of a second box, and still other plastic materials (e.g., 
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material indicated by circles 230 and 232) fall outside of both the first and second boxes. A 
different INBOX variable is associated with each of the three regions. By using these three 
variables, it is possible now to estimate more accurately the particular polymer property of interest 
regardless of which box or region the material being tested falls into. 

In one embodiment, if the sample is a polyethylene sample falling within the middle range, 
the MI regression model associated with that, range will automatically be selected by the 
corresponding variable model yielding the jargest value of the other two variable models, and that 
MI model is then used to achieve the most accurate MI measurement possible. Note that if the 
other INBOX variables yield values that are too large, no MI model will be selected and thus no 
MI estimation or prediction will be made in these cases. 

A plurality of regression models for a property other than MI (e.g., xylene solubles, 
density, rubber/oil content. FRR, or MF) can be developed during the calibration procedure (steps 
1 50-1 56 of FIG. 5A). For .example, two or. more regression models for rubber content (e.g., R, 
and.Rj) can be developed, b,this case, the discrete (e.g., two-valued or Boolean) variable 
regression model allows a prediction to be made about which of the two rubber content 
regression models should be used for any particular unknown sample. After the rubber content 
model is selected by using the discrete variable model, the selected rubber content regression 
model is applied to data derived from, the unknown sample to determine the rubber content of the 
unknown sample. 

Referring again to FIG. 5A. the steps performed generally are as follows. Samples with 
known types, properties, and quantities of target nuclei (including flow rates such as MI, FRR, 
and MF in plastics such as polypropylene, polyethylene, and polystyrene) are introduced into the 
NMR system (step 1 50). The FID curve is digitized via a flash converter of at least 12 bits 
accuracy and the result is stored in memory. The M-L iterative process is then used to derive the 
curve coefficients from the stored FIDsto a given Chi-squared error and iteration limit (step 152). 
The various non-linear transformations to be used are then determined to arrive at the various 
parameters (step 154). The types, properties, and quantities of target nuclei in the known samples 
are then related to the parameters by a regression against these transformed parameters to arrive 
at a regression equation or model for the discrete INBOX, variable or variables and for each 
physical property of interest including at least two separate models for at least one polymer 
property such as MI (step 1 56). These models are stored in the NMR system in memory. 
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Unknown samples can now be introduced into the NMR system. The FID curve for an 
unknown sample is recorded and digitized (step 158). The curve coefficients are then derived 
(step 160), and the parameters are calculated (step 162). Based only on the parameters of the 
unknown sample (and not on any external user-entered information), the INBOX regression 
model(s) is/are used to arrive at a particular value for the discrete INBOX variable(s) (step 164). 
This value (or values) corresponds to one of the polymer property (e.g.Jvfl) regression models, 
and that model is selected (step 164). In the case of more than two regression models for a 
particular polymer property of interest, the same number of discrete INBOX variable regression 
models are used, and each INBOX variable model corresponds to a different one of the property 
models. In the case of just two property models, only a single INBOX variable model is used, and 
each value of the two-valued .variable corresponds to a different one of the property models. In 
either case, the selected regression model for the property of interest is then applied to the 
parameters derived from the unknown sample to determine a value for that property of the 
unknown sample (step 166). The property can be estimated for other unknown samples by . 
returning to step 1 58 and introducing into the NMR system the next unknown sample, as 
indicated by a dotted line 168. 

Further details of the steps described generally above and with reference to FIG. 5 A are 
presented below with reference to FIGS. 6A and 6B. 

Note that it may be preferred not to use INBOX for some embodiments of the NMR 
system. . 

Referring to FIG. 6A, the first step 36 is to measure samples with known types, properties 
and, quantities of target nuclei, including flow rates (e.g!, MI, FRR, and MF) in plastics (e.g., 
polyethylene, polypropylene, and polystyrene). This data gathering (step 36) may be done on-line 
or off-line. The FID curve is digitized via a flash converter of at least 12 bits accuracy and stored 
in memory. The step 38 is to apply the M-L iterative process to derive curve coefficients from the 
FBDs to a given Chi-squared error. In step 40, the second order "x" variables (i.e., the 
explanatory variables) are formed. These x variables can include various parameters such as ratios 
of Y-axis intercepts, squares and cross products of these ratios, decay times, and temperatures. 
Higher order combinations of these parameters also may be calculated. These x variables .can be 
thought of as vectors in a multidimensional space where the space dimension is equal to the 
number of the explanatory variables. If there is no multicollinearity among the x variable, data, the 
vectors are orthogonal in this space (i.e., all dot products are zero). As multicollinearity 
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increases, the vectors move away from orthogonality. In the extreme case, there may be perfect 
correlation between two or more of the x variables and the vectors will lie on top of one another. 
An attempted regression analysis of such data would generate singular matrices. In step 42, x 
variable data with correlations above a certain threshold (e.g., 0.99) can be eliminated. 

The next step (step 44) is to choose a set of potential explanatory variables, the x data, 
from the M-L derived time equations including second and higher orders of these variables. They 
are chosen by a stepwise technique or some other .known technique. In a preferred embodiment, 
three different sets of x variables are selected and taken through the entire remaining steps (steps 
46, 48, and 50) and the set giving the best results is used. In a preferred embodiment, the best 
result is that which results in the jowest adjusted standard deviation of error on the degrees of 
freedom. One of the three different sets is composed of all the x variables. The second set is 
formed by the known stepwise technique of adding each new variable and determining if that 
variable helped and then continue adding those variables help. The technique is also applied 
in a backwise fashion where each previously added variable is retested in the presence of each new 
V ^ ,e \..^*^ rd .^ : i^^ ^.independent variables and variables with 
correlations between selected low and high limits, usually 0.2 to 0.96. 

The dependent, variables (the "y" variables) represent the property or characteristic of 
interest for the sample being tested. The y variables are related by a set of linear equations to the 
x or explanatory variables. 

The next step (step 46) is to perform a Principal Component Analysis (PCA). Every linear 
regression model can berestated in terms of a set of orthogonal explanatory (x) variables, where 
the new variables are linear combinations of the original x variables. The new x variables are 
called principal components and are orthogonal, thus eliminating the problem of multicollinearity. 
The regression model equation using the original explanatory variables is 

Y-Xp + u (Eq. 1) 

where Y is an n-by-1 column matrix of n observations, X is an n-by-p matrix of n observations on 
P explanatory variables, 0 is a p-by-1 column matrix of regression coefficients, and u is an n-by-I 
column matrix of residuals. If it is assumed that the expectation of u is 0, and that the expectation 
of uu',< wherein tf is the conjugate of u) equals the variance times the identity matrix, and that X 
and Y have been centered and scaled so that the XX' and YY' are matrices of correlation 
coefficients, then there exists a matrix C, satisfying 



WO 97/42516 



PCT/US97/05293 



-20- 

C(XX>)C=A (Eq. 2) 

and 

C'C=CC=I (Eq. 3) 

where A is a diagonal matrix with ordered Eigenvalues of XX* on the diagonal. The columns of 
C are the normalized Eigenvectors. 

A new set of explanatory variables Z may be formed by Z=XC. These are summarized as 

Y=X0+u=XCC f p+u=Za+u (Eq. 4) 
where the Z vectors are orthogonal. 

This process (step 48) of transforming the x data into z data produces a diagonal matrix A 
of Eigenvalues of the principal components. An Eigenvector of a squarematrix A of order n is a 
nonzero vector v where Av=Xv, and the scalar X is caJled an Eigenvalue. Eigenvalues may be 
calculated for matrix A from 

IA-A.IH) (Eq. 5) 

where I is the identity matrix, and the corresponding Eigenvectors v may then be found by solving 
(A-XI)v=0. The Eigenvalues are sorted numerically from the largest (top left of the diagonal) to 
the smallest (bottom right). If strongmulticollinearity exists, as it does for many of our 
explanatory variables, one or more of the lower right diagonal terms Will be very small compared 
with the others and these terms may approach zero. If the Eigenvalue is sufficiently close to zero, 
the value of the corresponding z transform of the x data is also essentially zero, as given by 

Z k =F k (xl,x2. . . ,xn)-0 (Eq. 6) 
where F k is the linear transform derived from PCA. 

The relationships of Eq. 6 are used to test each M-L curve fit to determine whether the x 
values obtained from M-L comport with those in the calibration set (derived from the known 
samples). For this test, a standard set of parameters are used. The set includes five variables: R ag 
(a ratio of Abragam amplitude/Gaussian amplitude); R« g (a ratio of exponential 
amplitude/Gaussian amplitude); t c (exponential T2); t, (Abragam T2); and t g (Gaussian T2), and 
their fifteen cross products for a total of twenty variables. The highest existing Z k (usually Z20) is 
chosen for the test, A seven sigma limit on the value of the selected Z k is used, and M-L solutions 
which result in Z k 's which are outside this range are rejected as M-L fitting failures. 
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The orthogonal explanatory (Z) variables are used in the regression to determine a model 
equation 48. Since the Z variables are orthogonal, the stepwise technique (or another known 
technique) is a reliable method for selecting terms for use in the regression model. 

The steps shown in FIG. 6A and described herein can be used to form a regression model 
for each property of interest. For example, a regression model for density, a regression model for 
xylene solubles, a regression model for MI, a regression model for FRR, and a regression model 
for MF can be formed via the steps detailed in FIG. 6A. Model making includes x and y outliers 
detection by F-statistic on x,y residuals and leverage calculation and model validation by 
calculating, for example, root mean squared error of cross-validation (RMSECV) according to 



RMSECV = 



i 



(Eq. 7) 



where I is the number of calibration samples, and y, is the input observation value of a polymer 
property for a sample i, and y, is the predicted value of the polymer property when a model is 
constructed without a sample i.. 

In order to improve the accuracy of the estimations of a particular property of interest 
(e.g., MI), two or more regression models for the same property can be formed. In some 
embodiments, a plurality of MI regression models.are formed, wherein each MI model is different. 
Also, a regression model for each discrete variable INBOX is formed. Each INBOX variable 
corresponds to a different one of the MI regression models. The accuracy of MI estimations is 
improved by forming and utilizing two different MI regression models, and INBOX is a two- 
valued or Boolean variable where one value corresponds to one of the MI regression models and 
the other value of INBOX corresponds to the other MI model. 

The regressionmodels formed by the steps of FIG. 6A are applied to unknown samples in 
the manner indicated by the flowchart of FIG. 6B Referring now to FIG. 6B, the FID of the 
unknown sample is measured (step 60), and an M-L analysis is performed (step 62) from which 
the Eigenvectors (i.e., Z k 's) are calculated (step 64). The selected Z k 's (i.e., those with 
sufficiently small Eigenvalues) are tested against a deviation limitation (step 66) where the 
detemunant of the vector divided by the standard deviation of that vector in the calibration data 
will cause a rejection of the M-L solution when the result is greater than a particular value. In 
some embodiments, that value is 7, and in other embodiments, the value is 5 or 9. 
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The Z k 's for those M-L solutions which pass the limit test (70) are regressed via the 
model regression equation to predict the value of the discrete variables INBOX (step 72). As 
described previously, this value corresponds to one of a plurality of regression models for a 
particular property such as MI. In some embodiments, INBOX is a Boolean variable, and two 
Ml regression models are provided. One value (e.g., 0) of INBOX corresponds to one of the Ml 
models, and the other value (e.g., 1) corresponds to the other MI model. As indicated at step 74, 
the estimated or predicted value of INBOX is used to select one of the plurality of regression 
models formed for one of the properties of interest (e.g., to select one of two MI models). With 
the appropriate MI regression model now identified for this unknown sample, that model is then 
used to predict an accurate value for MI (step 76). Step 76 involves performing the procedures 
indicated by steps 62, 64, and 66 in order to predict MI. 

Note that if the limit test is failed (68), M-L (step 62) is restarted with different starting 
assumptions and steps 64 and 66 are repeated. Also, should repeated failures occur with a given 
unknown sample, that sample is discarded, and a new unknown sample is taken and the steps of 
FIG. 6B are repeated. 

As stated hereinabove, a partial least square regression (PLSR) technique may be 
employed in the development of system models useful for predicting a value of a particular 
property of a polymer sample in real time Generally, the PLSR technique is helpful in improving 
calibration model accuracy by including x variations relevant for calibration of the polymer 
property of interest in early regression factors and ignoring later factors which may contain large 
but irrelevant x variations. 

The PLSR technique is basically a multivariate calibration method applicable to the x 
variable data generated in step 40 of FIG. 6 A. The PLSR technique uses the independent y 
variables actively during bilinear decomposition of the x data to generate a new set of fewer, y- 
relevant explanatory variables or "t factors". 

In an orthogonalized version of the PLSR technique, for each t, factor, where a - 1, 2, 
3...n, the following steps are performed on the centered x, y data, generic method steps being 
generally discussed in Multivariate Calibration by H. Martens and T. Naes (John Wiley & Sons, 
1991) at pp. 121 to 123. First, the loading weights, W„ are calculated using a least square 
regression technique and the local "model": 

x..i - y»-i W; + E (Eq. 8). 
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The vector W , is scaled to a length of one. Next, the factors, t„ are estimated using the local 
"model": 

x a ., - t a WS + E (Eq. 9). 
The NMR loadings, P„ are then calculated as follows: 

x,., = ( , />,» +E (Eq. 10). 
Next, the polymer loadings, q„ related to density, MI, etc. are estimated: 

y.-i= '\<7, + f (Eq. II). 

New x and y residuals are. then calculated as follows: 

x, = x 4 .,- -/, p,* (Eq. 12); 

. ,. . y» = y.-i - / . 9.. . (Eq. 13). 
The calculation steps of Eqs. 8-13 are repeated until all needed factors, t. , are calculated. 

Thereafter, vectors b „ and b are computed for "a" partial least squares factors to be used in the 
prediction, 

1 ; y= \b 0 + xb (Eq. 14), 

as follows: 

b.~ W (P'W)' 1 <j (Eq. 15); 
bo^y-x'b (Eq. 16). 

Lastly, unknown polymer parameters ys may be calculated as follows: 

yi = i 0 + Xi' b (Eq. 17), 

where x ; are the input data. 

Depicted in FIG. 8 is a flowchart of steps employing the PLSR technique developed 
calibration model. In step 3 1 0, FID curves for known samples are obtained, for example by 
recording, digitizing, and storing FID signals as discussed hereinabove. Next, the M-L iteration 
process is applied to the FID curves to derive component curve equations having a fit within 
desired Chi-Squared error, as depicted in step 320. Other iterative t chniques such as Gauss- 
Jordan, Newton-Raphson, or "steepest descent" may be employed as desired for a particular 
application; however, M-L is generally preferred. The M-L iteration process yields best fit 
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component curve equations corresponding to Abragam, Gaussian, and/or exponential curves as 
depicted schematically in FIG. 4. In step 330, relevant constants such as ratios, squares, and 
cross products of Y-axis intercepts, decay times and temperature from the derived component 
curve equations are calculated. Then, in step 340, the PLSR technique described with respect to 
Eqs. 8-16 is applied to develop a system calibration model. Model making includes x,y outliers 
detection by F-statistic on x,y residuals and leverage estimation, calculating RMSECV (Eq. 7) for 
model validation and determining the optimal number of PLSR factors (components). Other 
validation/optimization techniques may include F-statistic on regression and calculation of the sum 
of squares of observation residuals divided by the number of observations after "a" PLS 
components are calculated. Steps 350, 360, and 370 respectively entail generating an FID curve 
or relaxation signal for an unknown sample, deriving component equation curves therefrom, 
calculating respective equation constants, and applying the calibration model of step 340 to 
predict a value of the desired property of the unknown sample: 

While application of this particular PLSR technique has been demonstrated to be helpful in 
removing the impact of irrelevant x variations in calibration modelihg, the disclosure herein of the 
use of the PLSR technique is exemplary and not meant to be limiting in this regard. It is 
contemplated that other PLSR technique variations may be employed, such as an orthogonal 
loading technique, a Helland algorithm (LS. Helland, Communication in Statistics - Elements of 
Simulation and Computation, 17, 581-607 (1988)), or non-linear (including quadratic) versions, 
to develop calibration models of improved predictive accuracy and all such variations are 
considered to be within the scope of the appended claims. Further, multiple models according to 
this technique may be used in conjunction with the multiple model selection method described 
hereinabove. 

Referring to FIGS. 1 A, 1C, and 7, to establish effective industrial measurements with the 
NMR system , the following procedures are followed. A single FID curve is established to see if 
the sample area is clear (Quick FID) in an abbreviated cycle of attempting to establish an FID 
curve. If the sample region is not clear (N), measurement is interrupted to allow valve V2 to open 
and jets J and gravity to clear the region. A new Quick FID step establishes clearance. Then, a 
new sample is admitted by closing valve V2, opening valve VI, and making such adjustments of 
probe P and line LI as may be necessary to assure sample acquisition. Jets J adjust and stabilize 
the new sample. 
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An electronic signal processing apparatus baseline is established in 3 to 4 cycles, each 
having + and - sub-cycles with the addition of C+ and C- to detect a baseline offset and 
compensate for it. 

Further adjustment is established by coils 124 to adjust HO (i.e., resonance), and this is 
enabled by ten to twenty field check cycles of FID curve generation,. The C- FID is subtracted 
from the C+ FID (this process eliminates small baseline offsets) to obtain a workable digitized FID 
signal which has a maximum value at resonance. HO. is adjusted via coil current generator 1 17 and 
coils 1 24 until such maximum is achieved. In another embodiment, HO may then be changed to 
offset the system by a given amount of about 0. 1 to 3 KHz. 

Then one or more (usually five to one hundred) measurement cycles are conducted to 
obtain a useable measurement. Each of these five to one hundred cycles involves a modulated 
transmission/reception/flash A-D ; conversion and storage of data. The curves are then averaged 
for M-L curve fitting, and the above-listed.inteccepts and ratios are established. Similar cycles, 
often somewhat abbreviated, can be applied for Quick FID, field check, and baseline correction 
purposes. Each of the sub-cycles (i.e., + and,-) of each such cycle involves a capture and 
utilization of thousands of FID points in data reduction. Each sub-cycle occurs on the order of a 
second, and the number of such sub-cycles employed depends on the desired smoothing and 
signal-to-noise ratio (S/N). Generally, S/N improves in a square root relationship to the number 
of cycles. .... 

As noted in the above-listed related patents and applications, sample tube composition can 
distort readings. If glass is not used (and it is preferred to avoid glass in industrial usage), the 
replacement should not be a hydrocarbon plastic. However, fluorocarbons can be effective in 
several applications since signals from fluorine appear far from resonance. These signals can be 
distinguished from hydrogen at the levels of sensitivity required and if desired can be filtered or 
distinguished. In other cases of higher sensitivity measurements (e.g. , for gauging relative 
proportions of amorphous and crystalline species in mixtures thereof); the sample container should 
be glass or non-protonic ceramic. In some instances, however, fluorocarbon or reinforced 
fluorocarbon can be used acceptably for polymer measurements. In all such cases, the point is to 
avoid sample containers with species that can couple with transmitted energy and generate a FID 
decay curve mimicking the samples. 

While particular processing and modeling techniques are described herein, it should be 
realized by those of ordinary skill in the art that other techniques can be used instead of or in 
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addition to what is described. For example, neural networks and neural network techniques can 
be used with the NMR system to produce estimations of a variety of properties for a variety of 
unknown polymer samples. More specifically, a neural network may be employed to develop 
enhanced system models which utilize not only calculated FED component equation constants 
such as Y-axis intercepts, decay times, and temperature to predict sample properties, but also 
such process factors as product type, dwell time between manufacture and measurement, 
production rate, and reactor temperatures. 

A simple, exemplary neural network 400, also commonly referred to as a neural net, is 
depicted in FIG. 9. The net 400 consists of an input layer 402 having three input neurons I,, I 2 
and I 3 ; a single optional hidden layer 404 having three hidden neurons N,, N 2 and N 3 ; and an 
output layer 406 having a single output neuron O,. Any number of input, hidden, and output 
neurons I, N, and O may be employed as well as any number of hidden layers 404. Input to each 
hidden neuron N is a weighted sum of the input neurons I. ' For N, then, the input may be 
represented as follows: '••J-' " 

N, input = I, w,, + I, w 2 , + I, w 31 , . (Eq. 18); 

where w^ is a weighting factor specific to each pair of input and hidden neurons I, Nr The Ni 
input is then further processed by an activation function which may have any of a variety of 
mathematical forms. For example, a simple.aptivation function may be linear. A presently 
preferred activation function is of the form of a sigmoid function, mathematically represented as: 

N output = 1/(1 + e (Eq. 19). ' ' 

The sigmoid function includes both linear and non-linear regions, advantageously permitting the 
network 400 to accurately simulate very non-linear functions. Lastly, the output neuron Oi is 
calculated as the weighted sum of hidden neuron outputs as follows : 

Oi output = (N, output)(ow,) + (N 2 output)(ow 2 ) + 
(N 3 output)(ow 3 ) (Eq. 20); 

where ow NO is a weighting factor specific to each pair of hidden and output neurons N, O. 

In developing a system model according to the neural net technique, a plurality of data 
sets are empirically determined from known polymer samples. Each data set includes values for 
input neurons, I, corresponding to FID component equation constants (coefficients), measured 
process parameters, and/or other process variables. In another preferred embodiment; the input 



WO 97/42516 



PCT/US97/05293 



3 



-27- 

neurons, I, can be scores for optimal number oF PLSR or PC A factors, calculated from FID 
component equation constants in addition to measured process parameters and/or other process 
variables. Each data set also includes corresponding value(s) for output neuron(s), O, such as 
laboratory measured data of one or more properties of the polymer sample. For use in a 
representative application , input neurons I,, I 2 and I 3 may respectively correspond, for example, 
to a ratio of Y-axis intercepts, a decay time, and a.temperature. The resulting output neuron, O,, 
may correspond to the measured value of melt index or other property of interest. 

Upon receipt of the multiple data sets, the net 400 "learns" from the data provided using 
established correlative techniques such as back propagation. In simple terms, values for the 
weighting factors w^and own 0 are determined which best cause the output data (prediction), O, 
to match the corresponding laboratory-measured data of a polymer property (trial or training set). 
The minimization of the sum of squared errors between the prediction and the training set is 
accomplished using the steepest descent technique. Other minimization techniques such as the 
Gauss-Jordan, Newton-Raphson, Marquardt-Levenberg, or sequential quadratic programming 
(SQP) routine can be employed. 

Depicted in FIG, 10 is a flowchart of steps employing the neural network technique 
developed calibration 'model: In step 410, FlD curves for known samples are obtained, for 
example by recording; digitizing, and storing FID signals as discussed hereinabove. Next, the M- 
L iteration process is/applied to the FID curves to derive component curve equations having a fit 
Within desired Chi-Squared error; as depicted in step 420. Other iterative techniques such as 
Gauss-Jordan, Newton-Raphson, or "steepest descent" may be employed as desired for a 
particular application; however, M-L is generally preferred. The M-L iteration process yields best 
fit component curve equations corresponding to Abragam, Gaussian, and/or exponential curves as 
depicted schematically in FIG. 4. In step 430, relevant constants such as ratios of Y-axis 
intercepts (and possibly also squares and cross products of Y-axis intercepts), decay times, and 
temperature from the derived component curve equations are calculated. Then, in step 440, the 
neural network technique described with respect to Eqs. 1 8-20 is applied to develop a system 
calibration model, Model making includes validation by calculating RMSECV (Eq. 7). As 
discussed, relevant process data other than that developed from the derived equations may be 
employed advantageously to improve the correlative accuracy of the model. Steps 450, 460, and 
470 respectively entail generating an FID curve or relaxation signal for an unknown sample, 
deriving component equation curves therefrom, calculating respective equation constants, and 
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applying the calibration model of step 440 to this data and other process data to predict a value of 
the desired property of the unknown sample. 

In a particular example, a neural net was developed with fifteen input neurons, I, including 
neurons for values of derived component curve amplitude ratios, T2 time constants, ratios of time 
constants, etc. A single hidden layer was used with between five and fi^een hidden neurons, N, 
and a single output neuron, O. In some cases, for example where the power of a particular input 
neuron was a variable, the corresponding weighting' factor, "w Wl was substituted with products of 
input neuron values. In general, the neural net successfully learned from the data sets provided, 
being able to predict values of polymer properties of unknown samples not included in the 
learning data sets with accuracies corresponding to the previously described statistical methods. 
For example, unknown sample densities were predicted to within 0:0005 gms/cc of subsequently 
measured values and melt indexes were accurate to within about 7% to 10% of subsequently 
measured values. ; 

Not only are neural networks advantageous insofar as their ability to accurately predict 
unknown polymer sample properties of interest, less system computing power is needed for 
calibration since the data sets used for learning do not heed to reside in system memory at one 
time. Further, the neural network technique provides for substantial non-linear prediction, 
obviating squared and cross product terms in the input data set, further reducing the computing 
power required over other methods. Multiple models according to this technique may also be 
used in conjunction with the multiple model selection method described hereinabove. 

Whether a neural network, PLSR, or other technique is employed to generate a system 
calibration model, absent compensation for the phenomenon commonly referred to as u resin 
aging," predictive error may inherently exist. The amount of error depends on the parameter of 
interest as well as a variety of factors linked to the polymer manufacturing and measurement 
process. For example, in the case of a sample of polyethylene, initial sample crystallinity is 
affected by the manner in which the resin is cooled after formation. Depending on the 
temperature of the resin and the cooling n>te, the initial sample crystallinity may be close to or 
remote from an equilibrium crystallinity value for the material. The further from the equilibrium 
crystallinity value, the more the sample will tend to "age" over time as the resin relaxes and 
crystallinity approaches the equilibrium value. Further, the rate of aging is sensitive to resin 
temperature and time at temperature. A pelletized polyethylene sample, which is produced by 
heating and quench cooling the resin, will therefore tend to have an initial crystallinity which 
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deviates more significantly from the equilibrium crystallinity value than a counterpart powdered 
polyethylene sample. Absent the capability of the model to account for such variability, accuracy 
in predicted properties may suffer, sometimes dramatically. Development testing over an 
extended three hour time cycle to ascertain the magnitude of the problem has revealed variation in 
5 predicted melt index of a powdered polyethylene sample of up to 40.% to 50% and about an order 
of magnitude greater variability of a peptized polyethylene sample. While average time cycles 
prior to measurement routinely employed in measuring polymer properties of interest are 
significantly less and therefore predictive errors greatly reduced, these extended cycle values 
demonstrate the magnitude of the aging phenomenon and the need to account for such 
10 characteristics to achieve a more accurate predictive NMR measurement system. Accordingly, a 
resin age factor has been developed for use in the predictive model to reduce errors and enhance 
model accuracy in predicting the value of polymer properties of interest including crystallinity, 
density, and melt index, for example. 

Resin age has been shown to be a ^nction of at least one of the following three process 
5 parameters: elapsed cycle time ; between a predetermined initial time and the time of measurement; 
sample temperature at ; time pf measurement; and sample form By using these process parameters 
for known samples in- generating the system model, and similarly using readily measurable values 
for these process parameters for unknown samples in combination with the model, error in 
predicted values of melt index, for example, can be reduced from, for example, between about 
3 10% to 15% to about. 6% to 7%,, , . 

Shown in FIG. 11 is a flowchart ofsteps performed in accordance with an embodiment of 
the invention that uses a resin age parameter, either in the initial generation of the system model 
with known parameter value samples pr in the use of such a model with unknown parameter value 
samples or both. Reference may also be made to FIGS. 1 A-l C. 

Looking first to elapsed cycle time, any of several initial time points can be utilized. A 
• particularly convenient time for a preferred embodiment of an on-line NMR system is the time at 
which the sample is captured from the industrial process line IPL in step 500 and transferred to 
the heating chamber 170. A clock or counter, which may be integral with or external to the 
computer 106, is started in step 510 at time of transfer Clearly, an earlier clock initiation point 
may be utilized, for example at the time resin is first received in the storage bin from the reactor 
(not depicted). While applicable to an on-line NMR system as described herein, such an earlier 
initial time point may be particularly advantageous when using the NMR system in a non-real 
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time, off-line environment such as a laboratory where a substantial amount of time may elapse 
between initial transfer of the resin to the storage bin and ultimate measurement in the laboratory 
NMR system. 

As discussed previously with respect to FIGS. 1 A and IB, once in chamber 170, the 
sample is heated in step 520 by heating element 174 and temperature controlled gas stream 176 
until sample temperature reaches the desired mobility enhancing temperature, a process which 
may take several minutes. In step 530, inlet valve VI is opened and the sample is advanced into 
the critical sample region S2 where sample temperature is maintained by environmental control 
chamber 132 and air curtain containment jacket 137. Tuning of the magnetic assembly 1 16 to 
resonance with the sample in place follows in step 540. Once assembly 116 is tuned, the counter 
is stopped and elapsed time stored in steps 550 and 555, respectively. Thereafter, the excitation 
of coil 100 and excitation-precession of the proton content of the sample with subsequent 
relaxation/decay produces the FM signal resulting in the FID curve which is then measured in step 
560 and curve fit in step 570. Calculations are performed on the fitted component curves in step 
580 to derive the desired component curve equation constants used in combination with elapsed 
cycle time and other process data comprising the resin age factor in the model to predict the value 
of the parameter of interest. Lastly, the sample region S2 is emptied by air jets J in step 590 and 
the NMR system is ready to accept another sample. Alternate methods of sample transfer using 
the apparatus of any of FIGS. 2A-2D may readily be substituted for the transfer method described 
withrespecttoFIG.il. 

In practice, elapsed cycle time between initial transfer of the sample to the heating 
chamber 170 in step 500 and measurement in step 560 for an on-line NMR system may be on the 
order of about 60 to 380 seconds or more, during which time aging occurs. For example, heating 
of the sample to the desired mobility enhancing temperature may take from about 30 to 350 
seconds or more, depending on the initial temperature of the sample, and tuning of the magnetic 
assembly 1 16 may ordinarily take up to an additional 30 seconds or so. 

To further improve predictive accuracy, the resin age factor may also include sample 
temperature at the time of measurement, for example at about the step of modifying the 
precession. Between the time the sample is heated in step 520 and the time measurement occurs 
in step 560, some measurable amount of cooling can occur during magnetic assembly tuning with 
the sample in the critical region S2. The environmental control chamber i 32 and air curtain 
containment jacket 137 significantly reduce such cooling effects, thereby minimizing temperature 
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induced errors. However, when these integral thermal control systems are malfunctioning or 
inoperative, accurate predictions can still be made with a model employing the resin age factor 
based upon sample temperature. Alternatively, it is contemplated that less costly thermal control 
systems can be employed or the systems eliminated altogether if a sufficiently effective 
temperature compensation model is utilized. Sample temperature, fo.r purposes of calculation of 
the resin aging parameter, may be readily obtained from thermocouple 133 or infrared sensor 133- 
1, as depicted in FIG. 3. This value may be conveniently stored with the elapsed cycle time in the 
computer 106 or elsewhere as desired. 

By employing a resin age factor dependent on parameters such as cycle time, sample 
temperature, and/or sample form, a simplified system model can be developed and used to achieve 
improved accuracy predictive values of properties of interest as. described in the method steps of 
FIG. 12. Similar in methodology to the development of the system models depicted in FIGS. 8 
and 10, FID curves for samples with known, properties are obtained in step 610, for example, by 
recording, digitizing, and scoring FID signals as discussed hereinabove. Next, the M-L iteration 
process is applied to the FID curves to derive component curve equations having a fit within 
desired Chi-Squared error, as depicted in step 620. Other iterative techniques such as Gauss- 
Jordan, Newton-Raphson, or "steepest descent" : may be employed as desired for a particular 
application; however, M-L is generally preferred. The M-L iteration process yields best fit 
component curve equations corresponding to Abragam, Gaussian, and/or exponential curves as 
depicted schematically in FIG. 4. In step 630, relevant constants such as ratios, squares, and 
cross products of Y-axis intercepts, and decay times from the derived component curve equations 
are calculated. Next, process data such as elapsed cycle time from step 555 in FIG. 1 1, sample 
temperature stored at time of measurement, and sample form are acquired and the resin age factor 
calculated in step 640 and used with the equation constants from step 630 in step 650 to develop 
a system calibration model by any of the techniques discussed herein, including principal 
component regression, PLSR, and neural network techniques. Steps 660 through 690 
respectively entail generating an FID curve or relaxation signal for an unknown sample and 
deriving component equation curves therefrom; calculating respective equation constants; 
acquiring the resin age factor process data parameters and calculating the resin age factor 
therefrom; and applying the calibration model of step 650 to predict the value of the property of 
interest of the unknown sample. Clearly, the acquisition of process data and calculation of the 
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resin age factor may be accomplished at any convenient times and are not restricted to the specific 

order of FIG. 12, either in developing or using the system model. 

In general, use of the resin age factor formulated from process data reduces the number of 

discrete input terms required in the system model to achieve the aforementioned improved 
5 accuracy results. For example, the number of input terms for a system model for predicting melt 

index may be reduced from between about ten to fifteen terms down to about five to six terms. 

Similarly, the number of model input terms for predicting density may be reduced from about ten 

terms down to as low as three terms. Further, multiple modeis.utilizihg the resin age factor term 

may also be used in conjunction with the multiple model selection method described hereinabove. 
0 Variations, modifications, and other implementations of what is described herein will 

occur to those of ordinary skill in the art without departing from the. spirit and the scope of the 

invention as claimed. Accordingly, the invention is to be defined not by the preceding illustrative 

description but instead by the following claims. 

What is claimed is: >. 
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Claims 

1 1 A method for predicting a value of a property of a polymer sample in real time from an on- 

2 line production process, said method comprising: 

3 storing a model useful for predicting a value of a property of interest of a production 

4 polymer sample, said model generated by: 

5 a ) acquiring free induction decay curves for polymer samples having known values 

6 of said property of interest; 

7 b ) applying an iterative technique to derive respective component curve equations 

8 from each of said . free induction decay curves; 

9 c) calculating respective component curve equation constants; 

I 0 ■ . . °) acquiring process data for said known samples; and 

II e ) generating said model using said constants, said process data, and said known 

12 values; 

1 3 providing a production sample having an unknown value of said property of interest; 

1 4 applying a base magnetic field to the production sample to effect precession of nuclei of 

1 5 the production sample; 

16 modifying the precession; 

1 7 receiving a resulting relaxation signal representative of a free induction decay of nuclei of 

18 the production sample; 

' 9 digitizing said relaxation signal; 

20 applying an iterative technique to derive respective component curve equations from said 

21 relaxation signal; 

22 calculating respective component curve equation constants; 

23 acquiring process data for said production sample; and 

24 applying said constants and said process data to said model to. predict the value of the 

25 property of interest, wherein said process data comprises a resin age factor. 

1 2. The method of claim 1 wherein said resin age factor is a function of at least one parameter 

2 selected from the group consisting of cycle time, sample temperature, and sample form. 

13. The method of claim 2 wherein said cycle time comprises time lapse between about the 

2 step of providing a production sample and the step of modifying the precession. 
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4. The method of claim 2 wherein said sample temperature comprises temperature of a 
sample at about the step of modifying the precession. 

5. The method of claim 2 wherein said sample form is selected from the group consisting of 
powder form and pellet form. 

6. The method of claim 1 wherein said iterative technique applied to either of said free 
induction decay curve or said relaxation signal is selected from the group consisting of 
Marquardt-Levenberg, Gauss- Jordan, Newton-Raphson, and steepest descent. 

7. The method of claim I wherein said component curve equations are selected from the 
group consisting of Abragam, Gaussian, and exponential curve equations. 

8. The method of claim 1 wherein said component curve equation constants are selected 
from the group consisting of ratios, squares, cross products of Y-axis intercepts, and decay times. 

9 . The method of claim I wherein said property of interest is selected from the group 
consisting of xylene solubles, density, rubber/oit.content, melt index, flow rate ratio, melt flow, 
and crystallinity. 

10. The method of claim 1 wherein said production polymer sample is selected from the group 
consisting of rubber, polypropylene, polyethylene, and polystyrene. 

11. The method of claim 1 wherein said model generating substep (e) is accomplished by a 
technique selected from the group consisting of principal component regression, partial least 
square regression, and neural network techniques. . ' 

12. Apparatus for predicting a value of a property of a polymer sample in real time from an 
on-line production process, said apparatus comprising: 

means for storing a model useful for predicting a value of a property of interest of a 
production polymer sample; - - 

means for generating said model comprising means for: 

a) acquiring free induction decay curves for polymer samples having known values 
of said property of interest, 

b) applying an iterative technique to derive respective component curve equations 
from each of said free induction decay curves, 

c) calculating respective component curve equation constants, 
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1 1 d ) acquiring process data for said known samples, and 

1 2 e) 8 enerat ' n g said model using said constants, said process data, and said known 

1 3 values; 

14 means for providing a production sample having an unknown value of said property of 

1 5 interest; 

16 means for applying a base magnetic field to the production sample to effect precession of 

17 nuclei of the production sample; 

1 8 means for modifying the precession; 

19 means for receiving a resulting relaxation signal representative of a free induction decay of 

20 nuclei of the production sample; 

2 1 means for digitizing said relaxation signal; 

22 , m eans for applying an iterative technique to derive respective component curve equations 

23 from said relaxation signal; 

means for calculating respective component curve equation constants; 
means for acquiring process data for said production sample; and 
means for applying said constants and said process data to said model to predict the value 
of the property of interest, wherein said process data comprises a resin age factor. 

I 13. The apparatus of claim 12 wherein said resin age factor is a function of at least one 

parameter selected from the group, consisting of cycle time, sample temperature, and sample form. 
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The apparatus of claim 13 wherein said cycle time comprises time lapse between about 
when the production sample is provided and when the precession is modified. 

I 15. The apparatus of claim 13 wherein said sample temperature comprises temperature of a 
sample at about the time the precession is modified. 

1 16. The apparatus of claim 1 3 wherein said sample form is selected from the group consisting 

2 of powder form and pellet form. 

The apparatus of claim 12 wherein said iterative technique applied to either of said free 
induction decay curve or said relaxation signal is selected from the group consisting of 
Marquardt-Levenberg, Gauss-Jordan, Newton-Raphson, and steepest descent. 

18. The apparatus of claim 1 2 wherein said component curve equations are selected from the 
group consisting of Abragam, Gaussian, and exponential curve equations. 
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1 19. The apparatus of claim 1 2 wherein said component curve equation constants are selected 

2 from the group consisting of ratios, squares, cross products of Y-axis intercepts, and decay times. 

1 20. The apparatus of claim 1 2 wherein said property of interest is selected from the group 

2 consisting of xylene solubles, density, rubber/oil content, melt index, flow rate ratio, melt flow, 

3 and crystallinity. • * 

1 21. The apparatus of claim 1 2 wherein said production polymer sample is selected from the 

2 group consisting of rubber, polypropylene, polyethylene, and polystyrene. 

1 22. The apparatus of claim 1 2 wherein said means for generating said model utilizes a 

2 technique selected from the group consisting of principal component regression, partial least 

3 square regression, and neural network techniques. 
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